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Abstract 

We have developed a set of numerical tools for the quantitative analysis of defect dy- 
namics in quasiperiodic structures. We have applied these tools to study dislocation motion 
in the dynamical equation of Lifshitz and Petrich [Phys. Rev. Lett. 79 (1997) 1261] whose 
steady state solutions include a quasiperiodic structure with dodecagonal symmetry. Ar- 
bitrary dislocations, parameterized by the homotopy group of the D-torus, are injected as 
initial conditions and quantitatively followed as the equation evolves in real time. We show 
that for strong diffusion the results for dislocation climb velocity are similar for the do- 
decagonal and the hexagonal patterns, but that for weak diffusion the dodecagonal pattern 
exhibits a unique pinning of the dislocation, reflecting its quasiperiodic nature. 

1 Introduction 

It was realized long ago that dislocations play an important role in determining the mechanical 
properties of a crystal — whether periodic or not. As a consequence, ever since the discovery 
of quasicrystals much effort has been invested in the experimental and the theoretical study of 
their dislocations El) leading to interesting observations such as the special role played by 
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phasons in the motion of dislocations. 

Here we study the behavior of dislocations in a model system based on the dynamical 
equation of Lifshitz and Petrich [11 henceforth LP] , which is a modification of the well-known 
Swift-Hohenberg equation [H]. The dynamics of the LP equation is relaxational, dtp = —5J-/5p, 
driving a continuous 2-dimensional density field p{x, y, t) towards a minimum of the Lyapunov 
functional, or "effective free energy", 

Hp} = jdxdy {--/ + |[(v2 + i)(v2 + e)pf - \p^ + (1) 

yielding a dynamical equation of the form 

dtp = ep- c(v2 + 1)2 + q^fp + p2 _ p3_ (2) 

Among their results, LP showed that \iq = 2cos(7r/4), q = 2cos(7r/6), oi q = 2cos(7r/12) and e 
is sufficiently small, stable steady-state patterns are obtained (from random initial conditions) 
with long-range order and square, hexagonal, or dodecagonal symmetry, respectively. Although 
the LP equation cannot describe the dynamics of real solid-state quasicrystals and is probably 
more appropriate for soft matter or fluid phenomena such as Faraday waves, it offers a unique 
opportunity for the quantitative study of dislocations in a system — exhibiting both periodic 
and quasiperiodic long-range order — whose dynamics is exactly known. 

2 Injection and tracking of dislocations 

The Lyapunov functional in is clearly invariant under any translation or rotation of space. 
The steady-state solutions that are obtained are symmetry-broken ground states of J-^ whose 
Fourier transform has the form 

p(r) = j;p(k)e^^-, (3) 

where the (reciprocal) lattice L is a finitely generated Z-module, i.e. it can be expressed as the 
set of all integral linear combination of a finite number D of d-dimensional wave vectors. Li 
the special case where D, called the rank of the structure, is equal to the physical dimension d 
(here d = 2), the structure is periodic. 

2 



Any particular ground state p(r) of T is indistinguishable from a whole set of ground states 
that are related by so-called gauge transformations 

VkGL: /9'(k) = e2'^*^('^)p(k), (4) 

where x(k), called a gauge-function, has the property that x(ki -|-k2) = x0^i)+x0^2), possibly 
to within an additive integer, whenever ki and k2 are in L. As described in detail by Drager 
and Mermin gauge functions form a vector space V* of all real- valued linear functions on 
the lattice L. Because L has rank D, any linear function is completely specified by giving its 
values on D integrally-independent lattice vectors. The space V* is therefore a L>-dimensional 
vector space over the real numbers. The space V* contains, as a subset, all the integral-valued 
linear functions on L, denoted by L*. Gauge functions in L* leave the ground-state density 
invariant. Gauge functions that belong to the quotient space V* /L* take the ground state 
described by p into a different, yet indistinguishable, ground state described by some other 
density function p'. Thus, one can parameterize all the related ground states of .7-" on a simple 
D-torus — the order parameter space V*/L*. 

With this in mind we can easily construct an arbitrary dislocation in a 2-dimensional 
periodic or quasiperiodic density. As we traverse in a loop around the position of the dislocation, 
say the origin, we locally change the ground state p{r) by a gradually varying gauge function 
which winds around the i^^ direction of the D-torus Ui times, and returns back to the original 
point. This is most readily accomplished by using an angle-dependent local gauge function 
Xe(k) that assigns to the i^^ basis vector b^*) of L the value UiO. Thus, the most general 
dislocation is characterized by a set of D integers (ni . . . nr)), which for a periodic crystal 
reduces to the familiar d-dimensional Burgers vector. Not surprisingly, the set of all dislocations 
forms a rank-Z) Z-module — the so-called homotopy group of the D-torus 7 . Examples of 
dislocations — injected in this manner into the square and dodecagonal ground states of the LP 
equation — are shown in Fig. ^ after a short relaxation time. 

Once injected into the structure, the positions of dislocations are tracked numerically, as 
demonstrated in Fig. ^ by filtering individual pairs of Bragg peaks in the Fourier transform. 
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Figure 1: Top-left (in both figures): Snapshot of the numerical solution of the LP equation© 
showing a square and a dodecagonal pattern, a short time after a dislocation has been injected. 
Bottom-left: Fourier transform of the pattern. Note the fuzzy Fourier coefficients, containing 
the information on the angle dependent local gauge transformation. Bottom-right: A pair of 
filterred fuzzy Bragg peaks. Top-right: Inverse Fourier transform of the filterred peaks clearly 
showing the dislocation. 

and then performing an inverse Fourier transform to visualize the dislocations present in each 
individual density wave. This allows us to follow the positions of the dislocations in real time 
and obtain quantitative measurements of their velocities as described below. 

3 Dislocation dynamics under stress 

We apply external stress on a structure, containing a single dislocation, by squeezing it in a 
particular direction and then allowing it to evolve under the dynamics of the LP equation. We 
quantify the amount of stress by the change 6q in the wavenumber of the fundamental density 
wave in the direction of the applied stress, relative to its wavenumber in the relaxed steady- 
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Figure 2: Measurement of dislocation climb velocity under stress for large values of the diffusion 
constant c. (a) Position of a dislocation core in a dodecagonal pattern as a function of time 
under Lifshitz-Petrich (LP) dynamics; (b) Dependence of the dislocation climb velocity on 
applied stress in a stripe pattern, where v oc 5q^/'^ for Swift-Hohenberg dynamics and v oc Sq 
for LP dynamics; (c) Dependence of the dislocation climb velocity on the diffusion constant c 
under LP dynamics for different patterns. 
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state density. We have calculated elsewhere W, using the approach of Siggia and Zippelius [HI 
and of Tesauro and Cross 10 , that under such circumstances, with the parameters used in 
our simulations, the dislocation should climb with a velocity which is linear in the stress 6q 
and is proportional to the parameter c, playing the role of a generalized diffusion coefficient. 
This should be contrasted with the fact that for dynamics governed by the Swift-Hohenberg 
equation [31 the climb velocity is proportional to Sq^^"^, for the same simulation parameters. 

Fig. 121 summarizes our measurements of dislocation velocities, obtained using the automatic 
numerical procedure for tracking the dislocation, described above. Fig. El^a) shows the position 
of a dislocation as a function of time for different values of stress 5q, as measured under the 
dynamics of the LP equation for a dodecagonal pattern. For sufficiently large values of the 
diffusion constant c the dislocation climbs at a relatively constant speed, which indeed varies 
linearly with the applied stress, as shown in Fig. IHb). 

Fig- Etc) compares the measurements for different steady-state solutions of the LP equation. 
We find that the stripe pattern is the easiest for dislocation climb, with the climb velocity 
V ^ 4:5c6q, and that the square pattern is the most resistant, with v « 3c6q. Interestingly, 
the proportionality constants for the hexagonal and the dodecagonal patterns are very similar, 
with V ~ 25cSq. A possible explanation might be that the two triplets of wave vectors making 
up the dodecagonal pattern, act independently as two hexagonal patterns during the climb 
process. 

4 Dislocation dynamics under stress in the limit of weak diffu- 
sion 

As the value of the diffusion constant c decreases the local features of the pattern become 
important and the dislocation no longer climbs at a constant rate. In a periodic pattern, 
such as the square or hexagonal densities, the motion of the dislocation nearly comes to a 
stop at regularly spaced positions, or pinning sites, as shown at the bottom of Fig. |21 These 
sites correspond to positions where peaks in the density must first be annihilated before the 
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Figure 3: Bottom: Periodic pinning of a dislocation in a square and a hexagonal pattern for a 
small value of the diffusion constant c. Top: Irregular pinning of a dilocation in a dodecagonal 
pettern due to the lack of periodicity. 

dislocation can continue in its climb. In the quasiperiodic dodecagonal pattern a qualitatively 
different behavior occurs in which the dislocation is pinned at irregular intervals, as shown at 
the top of Fig. 1^1 It seems likely that because the quasiperiodic pattern contains different local 
environments the dislocation is more strongly pinned at certain sites than others, and these 
sites never quite repeat. This phenomenon has been observed in other models of dislocation 
motion by Mikulla et al. ^] and Pradkin [T^ . 
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